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Abstract. We study a recently observed gravitational microlensing peak in the V-band 
light curve of Q2237+0305A using a relatively simple, but highly consistent with the data 
(the best-fit reduced \ 2 is ver y close to 1), physical model. The source quasar is assumed 
to be a Newtonian geometrically-thin and optically-thick accretion disk. The disk has an 
arbitrary orientation, and both blackbody and greybody emission spectra are considered. 
When the electron-photon scattering plays a role, the greybody spectrum will be a simplified 
version of the exact one. In our model the microlensing variability result from the source 
crossing a caustic straight line. The main goal is to estimate the black hole mass and the 
mass accretion rate in QSO 2237+0305 as well as to discuss the power and the weakness of 
the technique, some possible improvements, and the future prospects from multifrequency 
monitoring of new microlensing peaks. We also put into perspective the new methodol- 
ogy and the results on the central engine in QSO 2237+0305. From the fitted microlensing 
parameters and reasonable dynamical/cosmological constraints, it is concluded that QSO 
2237+0305 harbours a central massive black hole: 10 7 M Q < M < 6 10 8 M©. While the 
information about the central dark mass is very interesting, the mass accretion rate is not 
so well constrained. The typical values of the disk luminosity/Eddington luminosity ratio 
are in the (1 — 20)e range, where e < 1 is the emissivity relative to a blackbody and the 
highest L/LEdd ratio corresponds to the largest deflector motion. Therefore, in order to 
verify L/LEdd < 1, a relatively small projected peculiar motion of the lens galaxy and a 
greybody emission seem to be favored. 

Ksy words. Gravitational lensing - Microlensing - Galaxies: nuclei - Quasars: general - 
Quasars: Q2237+0305 
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1. Introduction 



The optical continuum from a QSO could be originated from a geometrically-thin and optically- 
thick standard accretion disk. This standard scenario is a good candidate to explain most of the 
non-variable background component, whereas the fluctuations on different time scales may be 
caused by different mechanisms (accretion disk instabilities, supernova explosions in a circum- 
nuclear stellar region, and so on). The Newtonian model of a standard accretion disk around a 
black hole was introduced by Shakura & Sunyaev (1973). In this model, the released gravitational 
energy is emitted as a multitemperature blackbody radiation, where T s (r) is the temperature at 
radius r. Thus, in the absence of extinction, the emitted and observed intensities obey Planck laws 
B Vs [T s (r)] and B^[T(r)], respectively, where the emitted temperature and the observed one are 
related through the cosmological redshift of the source z s : T{r) = T s (r)/(1 + z s ). As the emit- 
ted temperature profile depends on the black hole mass and the mass accretion rate, the actually 
observed intensity profile I v (r) = e v B v [T(r)] is determined from the two physical parameters 
of the source. In a general situation we must take into account an extinction factor e v due to dust 
in the host galaxy of the quasar and any galaxy in the way to the observer including a possible 
lens galaxy and the Milky Way. A relativistic version of the Shakura & Sunyaev model was pre- 
sented by Novikov & Thorne (1973) and Page & Thorne (1974). The relativistic model should 
be a useful tool to describe the physics of the innermost layers of the disk. However, as we are 
interested in the emissivity of a wide optical source from an inner edge of several Schwarzschild 
radii up to an outer edge of 10 2 -10 3 Schwarzschild radii, in a first approach, we can ignore the 
relativistic effects and take the Newtonian picture. 

When a lensed quasar crosses a microcaustic, it is seen an important fluctuation in the flux 
of one of its images, i.e., a gravitational microlensing high-magnification event (HME) appears. 
The microlensing light curve of the involved image is basically given by convolving the observed 
intensity distribution with the corresponding magnification pattern. Therefore, the behaviour of 
the HME depends on the physical properties of the source quasar, and this fact makes way to a 
discussion on the structure of the source from the analysis of the prominent microlensing event 
(e.g., Rauch & Blandford 1991; Jaroszyriski, Wambsganss & Paczyriski 1992). In principle, if 
we assume a standard accretion disk as the source of the optical continuum, a measurement of 
the black hole mass and the mass accretion rate may be possible and, recently, several papers 
pointed that the study of an individual HME may be used to infer these parameters (Yonehara 
et al. 1998; Agol & Krolik 1999; Mineshige & Yonehara 1999; Yonehara et al. 1999; Shalyapin 
2001; Yonehara 2001). 

Very recently the GLITP (Gravitational Lenses International Time Project) collaboration has 
monitored the four images of QSO 2237+0305 in the V and R bands. In each optical band, 
the GLITP light curve for the brightest image, Q2237+0305A, traced the peak of one HME 
with an unprecedented quality (Alcalde et al. 2002). The global flat shape for the light curve 
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of the faintest image, Q2237+0305D, suggested that the intrinsic signal is globally stationary, 
and this result supported that the global variability in Q2237+0305A is exclusively caused by 
microlensing. The peak of the HME was fitted to the microlensing curves resulting from face-on 
circularly-symmetric sources crossing a single straight fold caustic, and it was obtained that the 
only source models totally consistent with the GLITP data are the Newtonian standard accretion 
disk and its simplified versions (Shalyapin et al. 2002). The uniform and Gaussian disks led to 
an excessively high value of the reduced chi-square, but the simplified variants of the Newtonian 
standard accretion disk gave excellent results: the best fits are in very good agreement with the 
observations and the upper limits on the source size are highly reasonable. To avoid a possible 
anomalous ratio of fitted background fluxes Fqr/Fqy, it could be needful to consider additional 
R light, which is emitted from an extended region (e.g., Jaroszynski, Wambsganss & Paczynski 
1992). This possible i?-band extended source is however irrelevant in the estimation of the com- 
pact source size as well as in the measurements of the central mass and the accretion rate from 
the technique presented in this paper. Unfortunately, from the framework used by Shalyapin et 
al. (2002), we cannot measure the physical parameters of the engine. One can only determine a 
relationship between the black hole mass (M), the mass accretion rate (dM/dt), and the quasar 
velocity perpendicular to the caustic line (V±). Even from some reasonable velocity range, there 
is no way to separately infer M and dM/dt. 

In this paper (Section 2), we introduce a novel expression of the microlensing light curve 
when a Newtonian standard source crosses a single straight fold caustic. This novel approach 
permits to break the degeneracy in the estimation of M and dM/dt, and so, to measure these 
parameters for a source with an arbitary orientation. In Section 3, from the U-band microlensing 
peak found by the GLITP collaboration, we obtain estimates of the black hole mass and the mass 
accretion rate in QSO 2237+0305. In Section 4 we present the main conclusions and put into 
perspective the technique and the new results on the central engine in QSO 2237+0305. 

2. Microlensing light curve for a Newtonian standard accretion disk near 
to a caustic straight line 

Yonehara et al. (1998) and Shalyapin et al. (2002) have previously studied the microlensing light 
curve associated with a Newtonian standard source that crosses a single straight fold caustic. 
Yonehara et al. (1998) presented the expected light curve when a face-on source (M = 10 8 M Q , 
dM/dt 0.13 M Q yr _1 ) is strongly magnified by a given caustic line, while Shalyapin et al. 
(2002) discussed the time evolution of the monochromatic flux for a face-on source (the inner and 
outer edges were assumed to be r in = and r out = oo, respectively) near to a generic caustic 
line. Using the theoretical light curve reported by Shalyapin et al. and an observed HME, one is 
able to obtain information on the size of the involved source quasar, provided that some interval 
of the quasar velocity V± might be inferred from observational data. However, their approach 
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does not permit a measurement of the two main parameters of the source: the central mass and 
the mass accretion rate. 

Here we like to show the behaviour of a HME caused by a generic Newtonian standard 
accretion disk (with finite inner and outer edges, and an arbitrary orientation) crossing a generic 
caustic straight line, and more importantly, a framework in which is possible to measure the 
physical properties of the black hole-accretion disk complex. We begin with a basic ingredient 
of the monochromatic radiation flux: the observed intensity profile. The observed intensity from 
a part of the disk at (r,ip) is given by 

/.(r) = e^ra] = ^ exp[w ; r(r)] _ i ; (1) 

where (Shakura & Sunyaev 1973) 

T(r) 



3 GM 



{l + z s )-\ (2) 



8-7T or 3 

The temperature profile depends on the black hole mass (M) and the mass accretion rate 
(dM/ dt), as well as the redshift of the source (z s ) and the inner edge of the disk (rj„). Two physi- 
cal constants are also involved: the Stefan constant (a) and the gravitation constant (G). As usual, 
ri n is assumed to be thrice the Schwarzschild radius of the black hole, i.e., r in — 6GM/ c 2 . Both 
the temperature and the intensity are circularly symmetric. Another basic ingredient is the magni- 
fication pattern. In our problem (magnification near to a straight fold caustic), the magnification 
law is simple and well-known. Taking a cartesian coordinate frame in which the caustic line is 
defined by the y-axis, the magnification of a piece of the disk at (x,y) will be (e.g., Schneider & 
Weiss 1987) 

A(x) =A Q + a c H(x)/Vx, (3) 

where H (x) = 1 at x > and H(x) = at x < 0. In the new coordinate frame, the magnification 
law is only function of the coordinate x, and we can easily relate this cartesian coordinate to the 
corresponding ones (r,tp) in the natural coordinate frame of the disk. If the centre of the source 
is placed at x = x c , we found a relationship between x and (r,ip): 

x = x c + r sin ip cos a sin j3 — r cos ip cos (3. (4) 

To derive Eq. (4) it was considered that the axis of the disk is inclined by an angle a in the 
plane of the observer's sky, and the node line and the x-axis are separated by an angle (3. The 
elemental flux from the surface element rdrdip can be now estimated in a direct way. In the 
absence of the lens, dfl* = D~ 2 rdrdip cos a is the solid angle subtended in the observer's sky 
by the piece of the disk, where D s is the angular diameter distance to the source. The true solid 
angle will be dil — A[x(r, ip; x c )]d£l*, and consequently, the elemental flux has an expression 
dF v (r,%j}]x c ) — D~ 2 I v (r)A[x(r,ij)] x c )]rdrdip cosa. Moreover, the integral of dF„(r,ip; x c ) 
over the whole disk gives the total monochromatic flux F„(x c ), i.e., 

F v {x c ) = K u -/ {1 + 1 }d*P, (5) 

J r . n exp[hv/kT(r)\ - 1 J A Q y> x (r,^x c ) 
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where K u = (2hv 3 /c 2 D 2 ) A a e u cos a, and r out is the outer edge of the disk, which (r out ) could 
be of the order of 100-300 inner radii (ji n ). 

In order to derive the microlensing light curve from Eq. (5), we firstly calculated the integral 
over ip. The angular integration led to 



f-2lT 



a c H(x) 



f [1 + ^^P]# = 2tt 
Jo A Q 
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1 + ^TT G ^ 



(6) 



where / = (cos 2 (3 + cos 2 a sin P) 1 ' 2 , q = x c /rf, and 



-l 

Inserting all these results into Eq. (5), one obtains 



a W =l ^p=. (7) 



F v {x c ) = 2nK v r U \l + -£°G(x c /rf)] J* ■ (8) 



In a second step, we used the new variable £ = r/ri„ as well as the trajectory of the centre of the 
source x c (t) = V±(t — to). As usual, Vj_ is the quasar velocity perpendicular to the caustic line, 
and t is the time of caustic crossing by the source centre. The final light curve has an expression 



(9) 



where 
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C=^=, (12) 



^IOSMqJ ^10 26 g s- 1 J • ( 3) 

In Eq. (9), ^ ma a; is the ratio between the outer radius and the inner radius, i.e., ^ max = r out /ri n . 
Moreover, in order to reduce the computing time, from some MATHEMATICA packages, we 
inferred an analytical approximation to the function G(q) (see Appendix A). The exact function 
is very well traced by the approximated one, with a typical accuracy of about 0.01 %. We note that 
the time evolution of the theoretical HME depends on the frequency, two chromatic amplitudes 
{A v , B v = BA V ), and four achromatic parameters (£ TOax , C, D, and to). 

In practice one must compare the theoretical law (9) with the brightness record of a QSO 
image during a HME. In the comparison, the frequency v will be the central frequency of the 
optical band in which the observations were made, and as mentioned here above, £, max = 100 
{Tout = 300 Schwarzschild radii) or £, max = 300 (r out ~ 1000 Schwarzschild radii) are plausible 
choices of the edge ratio. Therefore, at a given ^ max (e.g., 100), the goal is to estimate the values 
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of the parameters A v , £>„, C, D, and to by fitting the observed HME. Once the parameters C 
and D have been measured, from Eqs. (12-13) we make the relations: 

Hbh = 1(T 4 y ± (km s~ l )/[fC (day" 1 )], (14) 

Hac = H 2 bh/D\ (15) 

with (j, bh = M/1O 8 M and fi A c = (d M / dt) / 10 26 g s _1 . Although the factor related to the 
orientation of the source (/) cannot be directly obtained from the fit, its estimation is relatively 
simple. We simulated 10 4 pairs of values (a,f3) covering the whole range of possible orientations, 
and calculated the corresponding / values. From the simulations, we got and drew (Fig. 1) the 
probability distribution P(f). In Fig. 1 it is apparent that the bulk of the probability is associated 
with high values of /, and the maximum of the distribution is reached when / — > 1 . Only very 
particular orientations (nearly edge-on disks with a node line almost parallel to the caustic line) 
lead to small values of the factor /. From the probability distribution we inferred that < / > = 
0.84 and / = 0.84^o 02 at la confidence level (68.9 percent of the simulations; see the shading 
region in Fig. 1). Finally, in order to measure the central mass and the flow of matter, it is only 
required dynamical information. When that information is available, through Eqs. (14-15) one 
can estimate both physical quantities. 




Fig. 1. Probability distribution of the factor related to the orientation of the source. The shading 
region represents the 1 a confidence interval including both the mean value and the most probable 
one. 



Goicoechea et al.: The central engine in microlensed QSOs 7 

3. Physical properties of the accretion disk in QSO 2237+0305 

The GLITP collaboration has recently reported on the peak of a microlensing high-magnification 
event detected in both F-band and i?-band brightness records of Q2237+0305A, which was 
observed with an excellent sampling rate of about three times the OGLE one (Alcalde et al. 
2002). The event was discovered by the OGLE team (Wozniak et al. 2000), who have followed 
the evolution of the whole F-band HME. In a first paper, from the GLITP T^-band and i?-band 
peaks, Shalyapin et al. (2002) have analyzed the nature and size of the y-band and i?-band 
sources in QSO 2237+0305. To do this task, they used a family of source models. All the models 
led to theoretical microlensing curves with the same number of free parameters. This approach 
(comparison between models causing four-parametric theoretical curves) is very useful to discuss 
the feasibility of different scenarios. The authors concluded that only two rough versions of the 
Newtonian standard scenario are clearly consistent with the observed optical peaks. Thus we 
know that a Newtonian standard accretion disk is in agreement with the GLITP monitoring of 
Q2237+0305A, and in principle, it is possible to measure the central mass and the flow of matter 
(see Sect. 2). 

3. 1. Measurements of fC and D 

In our microlensing experiment, the functional relation presented in Eq. (9) is investigated by 
comparing it with the observed value of F v {t) at a given frequency and for various epochs 
t\, t2, £jv. A value of £ max is also assumed. The goal is to find the parameters A v , B v , C, 
D, and to of the theoretical microlensing curve which best describe the observational data. To 
make the fit, two different chi-square minimizations have been carried out. The most elementary 
procedure is to form a grid of points in the free parameters and evaluate the \ 2 function at each of 
these points. The point with the smallest value is then the minimum. Therefore, we have studied 
a reduced grid using some MATHEMATICA packages. The procedure is called RGMath and it 
is based on the ideas discussed in subsection 2.3 of Shalyapin et al. (2002). Basically, as we like 
to fit a theoretical curve that is linear in two parameters (A v and £>„), and non-linear in C, D, 
and to, the equations of minimization d\ 2 /dA L , = dx 2 /dB v = lead to analytic relations: 
A v = A V (C, D, to), B v — B V (C, D, to). In this scheme, the number of effective parameters is 
reduced to three, and one can work with a 3D grid instead of a 5D one. As the procedure is quite 
time consuming, we initially formed a 3D grid with moderate resolution, i.e., the size of the grid 
step was reasonable, but relatively large. So, from RGMath at moderate resolution, we inferred 
an initial best solution and some uncertainties associated with it. Another procedure is called 
DSFor and it is related to a version of the simplex method. We made a FORTRAN program to 
apply a downhill simplex methodology (Nelder & Mead 1965) to our problem. The simplexes 
will be the simplest geometrical figures in 5 dimensions having 6 vertices (points). The technique 
begins by choosing a point (vertice) in the 5D parameter space. Then the program generates an 
initial simplex, and after, repeated calculations are made while varying the vertices of the initial 
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Procedure 


^max 


A v (//Jy) B v (/Jy) 


fC (day- 1 ) 


D 


to (JD-2450000) 


RGMath 
DSFor 


100 
100 


2.78 1.37 
2.96 1.35 


2^+° 13 
u - zo -o.o6 

n 2i+ ' 14 

u -^°-0.08 


0.59 + _° f 3 
0.60l°;i 


1481.4 
1481.9 



Table 1. Parameters of the GLITP V-band microlensing peak. 

simplex in some way, until a local minimum is reached. The x 2 function may have different lo- 
cal minima around different points tested from the method, and consequently, to find the global 
minimum x 2 {min) we must make a 5D grid, and apply the downhill simplex technique to all the 
points in the grid. In the first stage with this alternative task (search of minima in the 5D param- 
eter space), the resolution was not very good. Once the initial best solution v* = (A*, B*, C*, 
D*, t$) is known (via RGMath or DSFor), it is possible to refine it and estimate accurate errors 
for the relevant parameters C and D. If the vector v = (A v , B v , C, D, to) of parameter values 
is perturbed away from v*, then x 2 increases. So we can draw high-resolution "parabolic" laws 
Xc (min) and x D [min), where, for example, Xc (min) represents the minima of the x 2 function 
for values of C around C*. The C range within which Xc(min) — x 2 (min) < 1 defines the la- 
confidence interval in the estimation of C, and the interval xb (min) — x 2 (min) < 1 is related 
to the lcr confidence interval in the estimation of D. The ka confidence interval in the generic 
parameter p corresponds to the region bounded by Xp(min) — x 2 (min) = k 2 (k= 1,2,...). 

The high-resolution fit to the GLITP y-band light curve for Q2237+0305A appear in Table 
1. The central frequency is of vy = 5.52 10 14 Hz, and we adopted £ max = 100. We note that 
the uncertainties presented in Table 1 are la intervals. Moreover, instead of direct measurements 
of the parameter C, la confidence intervals for the relevant factor fC (see Eq. 14) are quoted 
in Table 1. During an advanced stage of the project, we also obtained rough estimates of the 
errors in A v and B v . From DSFor we inferred A v - it\ uJy and B v - 1.5±^ /Jy. In Fig. 2 
we show the details to obtain the best values of C and the standard errors. The figure presents 
three Xc(min) trends, which were derived from RGMath at moderate resolution (dashed line), 
RGMath at high resolution (solid line), and DSFor at high resolution (dotted line). The la box 
corresponding to the RGMath (high resolution) procedure is plotted with a thin line. Looking 
both Table 1 and Fig. 2, two important conclusions appear: (1) the fit is stable against a change 
in the fitting method (RGMath or DSFor), and (2) a change from moderate to high resolution 
does not strongly perturb the parabolic law. Taking into account these last conclusions, we used 
the RGMath task (moderate resolution) to test the influence of the ( m „ value in the parameter 
estimation. We concluded that a change in £ max (from 100 to 300) does not modify the estimates 
of fC and D. In principle, a fit to the GLITP i?-band light curve for Q2237+0305A could be 
useful to improve the parameter estimation. However, the GLITP i?-band microlensing peak does 
not permit to infer the parameters C and D with errors similar to the previous ones (from the 
y-band peak). The ratio between the la i?-band interval and the lcr F-band interval is > 2.4 
for C and « 1.5 for D. Due to these relatively large errors from the i?-band brightness record, 
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Fig. 2. Behaviour of Xc(min) derived from the GLITP T^-band light curve for Q2237+0305A. 
The results with RGMath are shown by a dashed line (moderate resolution) and a solid line (high 
resolution). The dotted line was obtained with the DSFor task (high resolution). It is also depicted 
the la box corresponding to the RGMath procedure (high resolution). 

we only considered the estimates in Table 1. More properly, we took fC = 0.23+o'Jg day -1 
and D = 0.59^0 23 (RGMath solution) to determine the possible values of [i BIi and fi Ac . We 
remark that the accretion-disk model is consistent with the V"-band and i?-band light curves. The 
best-fit reduced x 2 is very close to one in both optical filters, and we can see the good agreement 
between best-fits and observational trends in Fig. 3. In the top panel, the observed V-band fluxes 
are compared with the best-fit from RGMath at high resolution (solid line), while in the bottom 
panel, the i?-band record is compared with the best-fit from RGMath at moderate resolution 
(dashed line). The fit in the R band led to best values of C = 0.375 day" 1 and D = 0.45, which 
are totally consistent with the blue parameters given in Table 1 . However, the i?-band data are 
relatively noisy in comparison to the ]/-band ones, and this fact does not permit a parameter 
estimation in the red band with uncertainties similar to the blue errors. 

3.2. Adoption of a dynamical range 

Kayser et al. (1986) showed that the effective transverse velocity of the source consists of three 
terms related to the transverse peculiar velocity of the source quasar ( v a ), the transverse peculiar 
velocity of the deflector ( Vd), and the transverse peculiar motion of the observer ( v ). The word 
"transverse" denotes "perpendicular to the line of sight", and the effective transverse velocity is 
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Fig. 3. Observed fluxes and best-fits. In the top panel (V-band), it is showed the best-fit from 
RGMath at high resolution (solid line). In the bottom panel (i?-band), it appears the best-fit from 
RGMath at moderate resolution (dashed line). 

given by 

v s v d D s v a D ds 

V = 1~> TZ TT + T, FT' W 

1 + z s 1 + z d D d 1 + z d D d 

where z d is the redshift of the deflector (lens galaxy), D d is the angular diameter distance to 
the deflector, and D ds is the angular diameter distance between the deflector and the source. In 
our case (QSO 2237+0305), the deflector is a nearby spiral galaxy (z d = 0.039 << 1) and the 
source is a far quasar (z s = 1.695). Due to these facts, one has that (1 + z s )/(\ + z d ) = 2.6 and 
D ds /D d ~ D s /D d <~ 10, and the particular configuration permits to neglect the term caused 
by the peculiar motion of the source. With respect to the motion of the observer, from the dipole 
observed in the cosmic microwave background (CMB), it is derived a peculiar velocity of the 
Sun: v obs = 369 km s" 1 in the direction (J2000) a = 1 l h l l m 57 s , 5 = - 7.°22 (Lineweaver et al. 
1996). This peculiar motion is almost perpendicular to the position for the system, being | v a \ w 
50 km s _1 . Because of the small observer's motion (as compared with the expected transverse 
peculiar velocity of the nearby spiral), we can also neglect the last term in Eq. (16). Therefore, 
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for QSO 2237+0305, the effective transverse velocity takes the very simple form 

V = -TT~W- < 17 > 

1 + z d D d 

We are indeed interested in the effective quasar velocity perpendicular to the microcaustic of 
interest. Thus, if i is the unit vector in the direction perpendicular to the caustic line, then 

V d D„ 

V± = V.i = 7-7^-7^, (18) 

1 + z d D d 

where v d = | v^.i\ is the amplitude of the projected peculiar motion of the deflector. 

Only one direct measurement of V = \ V\ is currently available. The measurement is an up- 
per limit on V, and consequently on V±, which was inferred from the analysis of the light curves 
of the system (Wyithe, Webster & Turner 1999). So we have decided to adopt a velocity range 
based on the relation (18) and the relevant observational data. Firstly, one must find plausible val- 
ues of Vd, and with this aim, we analyzed a catalog of galaxy data that is included in the electronic 
archives of the CDS (Centre de Donnees astronomiques de Strasbourg). The catalog is a part of 
the Mark III full catalog of redshifts and distances (Willick et al. 1997), and it contains redshifts 
in the CMB frame and both forward and inverse Tully-Fisher (TF) distances of 1355 spiral galax- 
ies (Mathewson, Ford & Buchhorn 1992). As it is well-known, the differences zcmb — ttf are 
the projected peculiar motions of the spirals, i.e., v g .n g , where n g are the unit vectors pointing 
towards the galaxies. In Fig. 4 we can see two very similar histograms showing two distributions 
of v g = | v g .n g \. The darkish line represents the results from the inverse TF distances {tjtf), 
while the other line traces the behaviour from the forward TF distances (tftf)- Inhomogeneous 
Malmquist bias-corrected distances tftf are quoted in the catalog, which were computed using 
the density field obtained through the IRAS 1.2 Jy redshift survey. From the distributions in Fig. 
4, we concluded that the mean value is of < v g > = 663 km s _1 . This average will be within 
the interval of allowed values of v d . As a lower limit we took 100 km s _1 , i.e., v d > 100 km 
s _1 , in agreement with the observed distributions [P(v g < 50 km s _1 ) < 6%] and discarding a 
conspiracy of nature (a situation in which both v Q = | v a .i\ and v d have a value < 50 km s^ 1 ). 
On the other hand, the distributions P(v g ) were obtained with a sample including field spirals 
and cluster spirals, and it is evident that amplitudes of the projected peculiar velocities in excess 
of 1000 km s _1 cannot be rejected. However, the lens galaxy is not placed inside a rich cluster 
or close to a group of clusters, and thus, one may assume an upper limit of 1000 km s _1 , i.e., 
v d < 1000 km s _1 . To derive a reasonable range for V±, the second basic ingredient is the factor 
D s / (1 + z d )D d . This factor depends on the redshifts (z d and z s ) and the matter/energy content 
of the universe (cosmological parameters). In this paper we consider the two standard flat cos- 
mologies: (1) O = 1, Ao = 0, and (2) f2 = 0.3, A = 0.7. As usual, f2o is the present density 
parameter and A is the present reduced cosmological constant. So, each plausible value of v d 
leads to two solutions for Vj_. Taking into account the whole range 100 < v d (km s _1 ) < 1000, 
the characteristic redshifts of the system, and the two standard cosmologies, it is found that 765 
< V± (kms- 1 ) < 10548. 
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Fig. 4. Probability distributions of the amplitudes of the projected peculiar motions of 1355 
nearby spiral galaxies (see Mathewson, Ford & Buchhorn 1992). The darkish line corresponds to 
the results from the inverse Tully-Fisher distances, and the other line represents the results from 
the forward Tully-Fisher distances, which were corrected for inhomogeneous Malmquist bias. 

3.3. Results 

Our results (la measurements of M and dMj dt for two cosmologies and different values of v d ) 
are presented in Table 2. We tested four representative values of v d (see subsection 3.2). The 
black hole mass is reasonably well constrained. For a universe with a zero cosmological constant 
and a very small projected motion of the deflector, the microlensing data suggest the existence of 
a central black hole having a mass of w 3.3 10 7 M . The most massive black hole is inferred for 
a universe with a nonzero cosmological constant (A = 0.7) and a very large projected velocity 
of the deflector. In that case, M w 4.6 10 8 M . Apart from these "individual" estimates, the 
main result is the global range for M. We concluded that the dark mass in the heart of the quasar 
should be larger than 10 7 M© and smaller than 6 10 s M Q . 

With regard to the mass accretion rate, the situation is quite intricate. Firstly, we cannot obtain 
individual lower limits, since all the la estimates are consistent with dM/dt = 0. Therefore, 
only typical values and upper limits are showed in Table 2. The fact that the la error bars on 
the accretion rate turn out to be compatible with zero is related to the (dis)abilities of the new 
technique. We fitted the parameters C and D (together with A u , B u , and in), and used Eqs. (14- 
15) to estimate the physical properties of the central engine: M and dM/dt. Considering the 
results in Table 1, fC is determined with a relative error of about 40 per cent, while the relative 



Goicoechea et al.: The central engine in microlensed QSOs 



13 



Cosmology v d (km s" 1 ) M (10 s M Q ) dM/dt (M© yr" 1 ) 
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3.041?;™ 
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1000 


4.59tkS 


277.7 « 734.4) 



Table 2. Black hole mass and mass accretion rate in QSO 2237+0305. 



error in D is of about 45 per cent. If we take a given value of the cosmological parameters and 
v d , then the relative error in M will be of about 40 per cent, i.e., similar to the relative error in 
fC (see Eq. 14 and the results on M in Table 2). However, unfortunately, Eq. (15) indicates us 
that dM/dt depends on both (/C) 2 and D 4 . Due to this fact, the accuracy in the determination 
of dMj dt is very poor, or in other words, the relative error in dMj dt is greater than 100 per cent. 

Secondly, a realistic scenario should incorporate the electron-scattering opacity, which could 
play an important role at inner regions of the accretion disk (Shakura & Sunyaev 1973). 
The local radiation energy flux is determined by the gravitational energy release Q gra v = 
(3/8Tr)(GM/r 3 )(dM/dt)[l — (rm/r) 1 ' 2 ]. In the outer regions, where the free-free processes 
give the main contribution to the opacity, it is formed a Planck spectrum. The radiation energy 
flux must be Q ra d = vT 4 , and from Q ra d = Qgrav we obtain the standard thermal law. In the 
inner regions of the disk, the electron-photon scattering presumably plays a role in the opacity, 
so these regions radiate less efficiently than a blackbody: I Vs \T s (r)] = ei n [v s , T s (r)]B Vs \T s (r)] 
with £i n [v s , T s (r)] < 1. If one defines a constant factor £j„ < 1 to be the emissivity relative 
to a blackbody, then the corresponding radiation energy flux will be Q ra d — £in&T 4 - At small 
radii the real disk could be hotter than the standard disk without electron-photon scattering, i.e., 
T s (r) = [(3/8tt) {GM/ar z )ell {dM/dt){\ - y/r in /r)] 1/4 . Therefore we can consider two sim- 
ple pictures: a global blackbody spectrum (e = 1 along all the disk), i.e., the standard one, and a 
global greybody spectrum (e < 1 along all the disk). This last simplification was suggested by 
Rauch & Blandford (1991), who introduced the greybody spectra in microlensing studies. When 
a generalized model is considered (e < 1), we must reinterpret the results on dM/dt as measure- 
ments of e~ 1 (dM/dt). In this way, the typical values of dM/dt vary in the interval (1 — 300)e 
Mq yr -1 . For e = 1, we derive a typical range of 1 - 300 M Q yr~\ while for s = 0.1, the typical 
mass accretion rates are 0.1 - 30 M & yr -1 . 
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4. Conclusions and discussion 



We showed the time evolution of a microlensing HME caused by a generic Newtonian 
geometrically-thin and optically-thick standard accretion disk crossing a generic caustic straight 
line. Given a lensed QSO, in order to estimate the mass of its central black hole and the ac- 
cretion rate, we can thus compare the theoretical law with a gravitational microlensing HME 
observed in some image of the far source. The technique was applied to the gravitational mirage 
QSO 2237+0305, whose brightest image experienced an important change in flux during the 
1999-2000 seasons (Wozniak et al. 2000). The peak of the dramatic variation in Q2237+0305A 
(Alcalde et al. 2002) led to very interesting information on the black hole in the core of the 
source and an interval of typical values for the mass accretion rate. The main source of uncer- 
tainty comes from our lack of knowledge about the peculiar motion of the lens galaxy in the 
direction perpendicular to the caustic line (v d ). To accurately determine the typical flow of mat- 
ter, it is also necessary to know the physical processes involved in the emissivity of the disk. 
Although the standard scenario does not incorporate the electron-photon scattering, in general, 
it could strongly distort the standard Planck spectrum at the innermost layers of the disk (e.g., 
Shakura & Sunyaev 1973; Malkan 1983). 

From la estimates of the relevant parameters in the theoretical microlensing curve, we in- 
ferred 1 a measurements of the black hole mass and the accretion rate for a set of Vd values, which 
represents the whole range of reasonable choices (see a detailed discussion in subsection 3.2). 
Taking into account all the individual estimates in Table 2, one obtains that QSO 2237+0305 
contains a massive black hole: 10 7 M < M < 6 10 8 M©. The information about the mass 
accretion rate is very much poor. The main result consists of a typical interval for dM/dt. This 
interval is of (1 — 300)e M Q yr^ 1 , where e = 1 if the free-free processes are dominant (standard 
source) and e < 1 when the scattering plays a role (e.g., Rauch & Blandford 1991). In spite of the 
fact that the range of typical determinations of dMj dt is not excessively broad, all the individual 
la measurements are in agreement with no accretion but a flow of matter of almost lOOOe M© 
yr _1 cannot be ruled out. On the other hand, one direct constraint on v d was recently reported by 
Wyithe, Webster & Turner (1999). This group (Wyithe, Webster, Turner and other colleagues) is 
involved in a project to interpret the QSO 2237+0305 microlensing light-curves, and as a part of 
the program, they derived an upper limit of v t < 500 km s _1 , where v t is the galactic transverse 
velocity (Q = 1, A = 0). Therefore, considering that v d < 500 km s _1 and v d > 100 km s _1 , 
we find a central value of v d = 300 km s _1 . For a flat universe with a zero cosmological constant 
and a projected motion of v d = 300 km s _1 , the typical source parameters are: M = 10 8 M© and 
dM/dt w 10eM© yr" 1 . 

Several modeling aspects can be improved, anyway. For example, one must include the ef- 
fects of general relativity. One may also incorporate a new theoretical law to the gallery of mod- 
els: the time evolution of a microlensing event generated by an accretion disk passing close to 
a cusp caustic. Efforts in both directions are now in progress, and we will consider these im- 
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provements in future studies of new QSO 2237+0305 microlensing peaks. Fortunately, using the 
present model and monitoring the events that occur in the components of the system, the fu- 
ture prospects are very promising. Thus, from 10 microlensing peaks monitored in ten different 
optical bands (blue bands are required to avoid the possible contamination by other sources of 
light as the extrapolated IR power law spectrum, and therefore, the possible perturbation of the 
chromatic amplitude A v which is defined in Eq. (10) for a blackbody emission; see, for example, 
Malkan 1983), we must be able to noticeably reduce the current uncertainty in the microlensing 
parameter D: from about 0.3 (see Table 1) to 0.03. In a similar way, if each peak is recorded 
at 10 frequencies, the error in each C value will be lowered in a factor of about 3, i.e., from a 
mean uncertainty of w 0.1 day -1 to w 0.03 day -1 . As a result of new detailed monitoring pro- 
grams, accuracies of 5-10% in the microlensing parameters C and D can be easily achieved in 
the next years. If a global dataset is made for each microlensing peak (tracing the time evolution 
of the spectrum), and it is fitted to a model with six free parameters: A, B, C, D, to, and 7 ( A v 
= Av z+1 , B v = Bv s+1 ), the accuracy could be even better than a few percent. However, as it 
was previously remarked, we need a robust estimation of the lens galaxy motion to accurately 
measure the black hole mass. This "dynamical problem" is the only pitfall to obtain a robust 
measurement of the amount of dark mass in the centre of QSO 2237+0305. From another point 
of view, if an independent estimate of M were available, we would deduce the effective quasar 
velocity V± involved in the different microlensing events as well as the value of e^ 1 (dM/dt). 
Through a reasonable collection of V± data, one may derive the effective transverse velocity of 
the quasar (V), which is a basic piece in some microlensing analyses. Up to now we focussed 
on the parameters (C, D) and the properties of the source, however, the parameters (A v , B v ) 
are also relevant for other studies. For example, if it is available a multiband monitoring of a 
microlensing peak, we may estimate the ratios A v u\/A v tj\. These ratios will be directly related 
to the extinction ratios e^/e^j), and consequently, to the global (host galaxy + lens galaxy + 
Milky Way) extinction law. 

Finally, we wish to put into perspective the new technique to determine the physical parame- 
ters of a far quasar and the measurements for QSO 2237+0305. Evidence for central black holes 
in local galaxies comes from small rotating gaseous and stellar disks discovered with the Hubble 
Space Telescope, VLBI observations of H2O masers in Keplerian rotation in some galactic nu- 
clei, and near-infrared measurements of the radial and proper motions of stars in the cluster at 
the centre of the Galaxy (e.g., Ford et al. 1998). Analyzing these gravitational signatures of mas- 
sive black holes, the typical values of the black hole masses vary between ~ 3 10 s M Q for the 
Milky Way (Genzel et al. 1997; Ghez et al. 1998) and - 3 10 9 M Q for M87 (Ford et al. 1996; 
Macchetto et al. 1997). The central dark mass can be also measured in far AGNs using two clas- 
sical methods: analyses of the blue and ultraviolet regions of the rest-frame spectra (e.g., Malkan 
1983; Sun & Malkan 1989) and observations of the time delays between the variations of the 
continua and the variations of the broad emission lines (e.g., Wandel, Peterson & Malkan 1999). 
The reverberation technique provides an estimate of the size of the broad-line region (BLR) as 
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well as the central mass of a quasar at redshift zqso- If the line-emitting gas is gravitationally 
bound to the central black hole, the virial theorem implies that M oc Rblr& 2 , where a is the 
FWHM of the emission profile of the BLR at distance Rblr- Therefore, measuring the lag r 
between an intrinsic event at the central continuum source and the reply (reverberation) at the 
BLR, one obtains the size Rblr = ct/ (1 + zqso) an d the dark mass M (e.g., Vestergaard 2002 
and references therein). Apart from the reverberation mapping technique, there is a straightfor- 
ward method to get the parameters M and dM/dt. A comparison between the blue-ultraviolet 
excess flux of a QSO (over the extrapolated infrared power-law spectrum) and the theoretical 
integrated disk spectrum permits to deduce the two physical quantities (e.g., Malkan 1983). In 
a pioneer work, Shields (1978) found M = 10 9 M Q and dM/dt = 3 M Q yr _1 for a Newtonian 
blackbody disk modelling the central engine in 3C 273. To make the fits, Malkan (1983) used 
a relativistic standard accretion disk. The theoretical model included the effects of general rel- 
ativity, but the inclination and scattering atmosphere were ignored. For a sample of 6 quasars 
and assuming nonrotating black holes, the "big blue bumps" can be fitted by spectra of disks 
around massive (~ 10 8 Mq) black holes. The author also concluded that the six high-luminosity 
quasars are emitting energy at approximately their Eddington limits. The effects of disk inclina- 
tion have been more recently considered by Sun & Malkan (1989), while the effects of electron 
scattering were discussed by Wandel & Petrosian (1988). In this paper we applied an alternative 
methodology to determine (M, dM/dt): the monitoring of a lensed quasar when a microlensing 
HME occurs (e.g., Yonehara et al. 1998; Agol & Krolik 1999). The method can be seen as a 
complementary tool to prove the existence of massive dark objects and accretion disks in the 
centre of far QSOs. Taking as reference values M = 10 8 Mq and dM/dt lOe Mq yr _1 (see 
here above), the estimate of the central dark mass in QSO 2237+0305 is consistent with data of 
other galaxy nuclei. However, the accretion rate is relatively large for a blackbody disk (e = 1). 
As L/LEdd (dM/dt) /M, where L is the total disk luminosity and Lsdd is the Eddington 
luminosity, one obtains L/LEdd ~ 3 for e = 1, L/LEdd ~ 1 for e = 0.3, and L/LEdd ~ 0.3 for e 
= 0.1. Therefore, to avoid meaningless results (L > L Edd), the electron-photon scattering must 
play a role. 
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Appendix A: approximation to the function G(q) 

The function G(q) is denned from Eq. (7) of the main text. For q > 1, we can rewrite G(q) as 




(A.l) 
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and expand the factor (1 + x) -1 / 2 , x = —y/q, — 1 < x < 1. Then we derive an approximated 
law 

n , tt(1155 + 1680g 2 + 3072g 4 + 16384 g 6 ) 

G(q) ~ 16384gi3/2 ' (A ' 2) 

which works very well at q > 2 (relative deviations less than 0.025%). However, the previous 
expansion does not work so well in the interval 1 < q < 2. Taking u = 1/ \Jq — y and dv = 
dyj \J\ — y 2 , G(q) will be equal to 

g <"»^(75 L t + v5Tt) +A ' ( ' ) ' <AJ) 

with 

y _ i . 1 f +1 , arcsiny , . JN 

j/.,*!^- (A - 4) 

Expanding the function arcsiny, we can approximate K{q) (and thus, G{q)) as 

At q ~ 1 the first approach (based on the expansion of (1 — y/q)" 1 / 2 ) underestimates the true 
behaviour, while the new scheme (based on the expansion of arcsiny) overestimates G{q). Due 
to this fact, at 1 < q < 2, G(q) is fitted to a law incorporating both trends. The final result is 

where a = 2.89283135146944, oi = 0, a 2 = 125.8537652305732, b = 
-102.1055508185553, b l = 227.8006617807037, b 2 = -125.694998412877, c = 
-202.3343397590727, a = -93.3006092335309, and c 2 = 124.1342775121369. The approx- 
imation (A. 6) is very accurate (relative deviations < 0.01%), even at q values very close to the 
singularity (q = 1). 

For q < 1, it is more difficult to obtain useful approaches to G(q). After some tests, we 
considered the combination of a smooth law (i.e., d n + diq + d 2 q 2 + ...) and a trend that is 
singular at q = 1 (i.e., (e + e\q + e 2 q 2 + — q). At — 1 < q < 1, from the combined law 

we infer a good fit 

G(g)^Y.drf + Y ^=» eiq \ (A.7) 

where d = 353.7491560566996, d 1 = -29.22210584958317, d 2 = -285.8452207779682, 
d 3 = -296.6766633600129, d 4 = -77.10649562470328, d 5 = 236.0210288747884, 
d 6 = 329.3059620964406, d 7 = -240.8295381320522, d 8 = 18.61151978935351, 
e = -351.1273028283772, e x = 205.3847185910097, e 2 = 315.1612274462573, 
e 3 = 171.9981477672829, e 4 = -95.1705826992738, e 5 = -320.8255310602542, e 6 = 
-243.8876547053204,e 7 = 415.2244998079605, and e 8 = -96.7574150920885. The relative 
deviations are again less than 0.02%. Of course G(q) = at q < —1. 



